Anderson localization of a Bose-Einstein condensate in a 3D random potential 
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We study the effect of Anderson localization on the expansion of a Bose-Einstein condensate, 
released from a harmonic trap, in a 3D random potential. We use scaling arguments and the self- 
consistent theory of localization to show that the long-time behavior of the condensate density is 
controlled by a single parameter equal to the ratio of the mobility edge and the chemical potential 
of the condensate. We find that the two critical exponents of the localization transition determine 
the evolution of the condensate density in time and space. 
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Anderson localization |l[ is an ubiquitous phenomenon 
which occurs in the propagation of waves (electrons in 
solids microwaves @, light 0], ultrasound in 
random media. Multiple scattering from random obsta- 
cles and the resulting destructive interference suppress 
wave propagation and can lead to exponentially localized 
cigenstates. 

Recently there has been much interest in the possibility 
of observing Anderson localization of Bose-Einstein con- 
densates obtained by trapping and cooling bosonic atoms 
[S IE 0, Hi • A Bose-Einstein condensate is characterized 
by a macroscopic occupation of a single quantum state Q 
and hence exhibits quantum, wave-like behavior despite 
its macroscopic size. Atomic Bose-Einstein condensates 
subjected to random external (optical) potentials are po- 
tentially good candidates for observing Anderson local- 
ization of matter waves. Up to now, the experiments have 
focused on ID configurations where all single-particle 
eigenstates are localized. In a typical experiment, the 
condensate is created in an optical or magneto-optical 
trap. The trap is then turned off and the condensate is 
allowed to expand. 

In this Letter we study the expansion of the Bose- 
Einstein condensate in a 3D random potential. Unlike in 
ID, a critical energy (the mobility edge e c ) exists in 3D 
which separates extended and localized states. An eigen- 
state is extended (localized) if the corresponding energy 
is larger (smaller) than e c . When the condensate is re- 
leased from the trap, the atoms achieve kinetic energies 
up to the chemical potential fi of the trapped conden- 
sate. For weak disorder e c < /i, and a fraction of atoms 
diffuses away, whereas the remainder is localized, as was 
pointed out in We show here that, surprisingly, even 
for strong disorder e c > [i only a fraction of the conden- 
sate will be localized. We study the full dynamics of the 
condensate expansion by accounting for weak localization 
at energies e > e c , strong localization at e < e c , and crit- 
ical behavior around the mobility edge. Our main result 
is that the effect of disorder on the expansion of the con- 
densate is controlled by a single parameter e c //i, and that 



Anderson localization plays an important role even when 
the chemical potential of the condensate fj, is much larger 
than the mobility edge e c . We show that the behavior of 
the average condensate density n(r, t) at large distances 
r and long times t is governed by the critical exponents 
v and s of the localization transition. This could provide 
a direct way to measure these exponents. The density 
of the localized part of the condensate n(r, oo) does not 
decay exponentially with r, as one could have expected, 
but follows a power law. 

Consider a Bose-Einstein condensate oi N 1 atoms 
of mass m trapped in a 3D spherically-symmetric har- 
monic potential VL,(r), characterized by the trap fre- 
quency u, to which we add a Gaussian uncor related ran- 
dom potential V(r): B(Ar) = V(r)V(r + Ar) = uS(Ar), 
where the horizontal bar denotes averaging over an en- 
semble of realizations of the random potential V(r). The 
wave function ip(r, t) of the condensate obeys the (mean- 
field) Gross-Pitaevskii equation: 
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where g measures the strength of repulsive interactions 
between atoms. We consider a situation in which the ran- 
dom potential V(r) is initially absent [EQ- A sudden 
turn off of the confining potential 14; (r) then gives rise 
to a rapid expansion of the condensate. During a time 
t ^ to = 1/lo the initial potential energy of the trapped 
condensate is converted into kinetic energy [E HE] ■ The 
kinetic energy per atom becomes typically of the order 
of the chemical potential /z of the trapped condensate. 
The calculation that we present below applies both in 
the regime of strong (fx S> foui) and weak (// ~ ftio) inter- 
actions. Whereas the former is realized in current exper- 
iments Q, the latter can, in principle, be reached using 
magnetic Feshbach resonances (llj . 

The random potential is switched on at a time t > to- 
This creates a new energy scale at the mobility edge e c . 
At this stage the kinetic energy of the condensate be- 
comes much larger than its interaction energy, and we can 
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set g — in Eq. JT]). The neglect of the nonlinear term 
in Eq. ([1]) for times t > to can be justified in the absence 
of disorder [§] and has been validated in the presence of 
a random potential by numerical simulations in ID 
Based on this, we expect that the neglect of interactions 
at long times t > t should be valid in 3D as well, at least 
when the interactions are not too strong. The analysis 
of the complete nonlinear equation ([T]) is a formidable 
task that falls far beyond the scope of this Letter. Even 
in ID controversial numerical results have been reported 
recently [e| • We emphasize that even though we neglect 
interactions between atoms for t > to, we fully take them 
into account at earlier times t < to and the strength of 
interactions g enters our final results through fi. 

The experimentally relevant quantity is the time- 
dependent condensate density n(r,t) = |^(r,t)| 2 . For 
large distances r 3> t (where £ = h A/ n/um 2 is the mean 
free path [l]|) and at long times t ^> H/e c , we derive the 
following expression for the condensate density, averaged 
over an ensemble of random realizations of V(r), 

j^f\m\ 2 J_^deA(k,e)P e (v,t). (2) 

Here P e (r,t) is the probability density to find a particle 
of energy e, initially located at the origin, in the vicinity 
of r after a time t ("probability of quantum diffusion" 
[3]), A(k,e) is the spectral function, and 0(k) is the 
Fourier transform of the wave function <fi(r) of the con- 
densate after the first stage of expansion (i.e. after a 
time t > t ). Physically, d 3 k|0(k)| 2 /(2-7r) 3 is the num- 
ber of atoms with a momentum around Hk at the time 
when the random potential is switched on, and deA(k 7 e) 
is the probability to find an atom with energy around e 
among all atoms with momentum hk, in the presence of 
the random potential. In terms of the Green's function 
of the Schrodinger equation G(r,t), the Fourier trans- 
form of P e (r, t) with respect to time t is P e (r, fi) = 
G(r, e + fi/2)G*(r, e - fi/2)/27ri> e , where v t is the den- 
sity of states, and A(k, e) = — ImG(k, e)/ir. In previous 
work Q, a free-space expression A(k, e) oc 5(e — ek) was 
adopted with et = h 2 k 2 /2m. This appears to be a bad 
approximation for energies near and below the mobility 
edge, where the uncertainty in energy is large due to dis- 
order. 

To find the spectral function A(k, e) we go beyond the 
first-order Born approximation [fj, LL4|, and use the so- 
called self-consistent Born approximation. We solve self- 
consistcntly the equations for the Dyson Green's func- 
tion G: G(k,e) = [e — ek — S(k, e)]" 1 , and the self-energy 
S: S(k,e) = / d 3 k'B(k - k')G(k', e). For the uncor- 
rected random potential we obtain G(k, e) = [e — e* — 
e c /4 — £k + «S/2r e ] _1 for e > e*. The mean free time is 
r t = nh 4 1 ' \f2rr?l 2 u\J(. — e*. The edge of the spectrum e* 
depends on a cutoff needed to regularize the divergence 
of the integral over k' in the equation for E. The value 
of e* is not important for the rest of our analysis because 



all relevant quantities depend on e — e*. The mobility 
edge e c is assumed to obey the condition k(e c )£ = 1 and 
is located at e c = e* + h 2 /2m£ 2 . 

Another important ingredient of Eq. @ is the proba- 
bility of quantum diffusion P € (r, t). At large r»f and 
fi <C e its Fourier transform can be found in the hydro- 
dynamic limit of quantum transport theory: P e (r, fi) = 
exp[-r^-iCl/D e (Cl)}/4TTD e (Cl)r 0, where D e (Cl) is the 
dynamic diffusion coefficient. As we will see below, the 
large-r behavior of n(v, t) is dominated by atoms with 
energies e > e c , which justifies the use of the above ex- 
pression for P e (r,Cl) for t ^> H/e c . 

Finally, the momentum distribution of the expanded 
condensate at the time when the random potential is 
switched on is assumed to be given by |(/>(k)| 2 oc l — k 2 /k 2 
for k < fc M = y / 2mfi/h and otherwise. This expres- 
sion follows from the dynamic scaling for t ^ to 
and /i > Sw. In the absence of interactions fj, ~ Hu> 
is a small energy scale. If < e c , Eq. ([2|) reduces to 
n(r,t) = N J^ ^ deA(0, e)P e (r, t) and the precise profile 
of |0(k)| 2 has no importance. 

In the following we will discriminate between localized 
(e < e c ), diffusing (e > e c ), and anomalously diffusing 
(e ~ e c ) atoms. In the limit of very long times atoms with 
energies above the mobility edge have diffused away, and 
only localized atoms contribute to the condensate density 
n(r, t — > oo) at any finite distance r. For these atoms we 
can set D t {Cl) = -iCl£ 2 (e) [H, [H] . Inspection of Eq. © 
then reveals that n(r, oo) is determined by the critical 
dependence of the localization length £ on the energy e 
in the vicinity of the mobility edge e c : £(e) oc |e — e c | ~ v . 
We obtain 

N ( £\ ljv 
n(r,oo) oc /(e c //x) — x - , (3) 
r 3 \r ) 

where f(x) oc x 3 / 2 for x <C 1 and f(x) ~ const for x 3> 1. 
This important result demonstrates that the critical ex- 
ponent v can be determined from the spatial profile of 
the condensate density at long times. While the value 
of v is not known for continuous disordered potentials, 
numerical solutions of the Anderson tight-binding model 
yield v ~ 1.5 [l8j|, the self-consistent theory of localiza- 
tion predicts v = 1 [l9| ]. and a claim of v ~ 0.5 has been 
recently made for light in Ti02 powders [20I ]. 

The stationary density profile ^ may take a long 
time to be reached, especially for large distances r. It 
is therefore important to look for dynamical signatures 
of Anderson localization. The complete density of the 
condensate can be represented as a sum of the station- 
ary part, considered above, and a time-dependent part: 
n(r,t) = n(r, oo) + 8n(r,t). At a given r 3> I and for 
short times, 5n(r, t) is dominated by the fastest atoms of 
the condensate. For weak disorder (e c -c m)j these atoms 
typically have kinetic energies of order fi, well above the 
mobility edge. Hence, they are almost unaffected by lo- 
calization effects and diffuse with the "classical" diffusion 
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coefficient D e (fl) ~ = £ 2 /3r e . The integral over en- 
ergies in Eq. (|2|) can be then evaluated using the saddle 
point method. For a given distance r, Eq. ([2]) reaches a 
maximum at the "arrival time" 
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This result is equal to that found for a quasi- 
monochromatic wave packet with central energy fi. 
At long times and for any e c //i, the dynamic part of the 

condensate density is dominated by the critical behavior 

> 
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of D e (fl) at energies e > e c : D c (fl) oc |e — e c 
The contribution of the critically diffusing atoms to the 
atomic density is 



N 

Sn(r,t) oc /(e c //iW 



l/s 
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(5) 



where the function f(x) is the same as in Eq. © and 
where s > 2/3 is required to assure convergence of the 
integration over energies. Equation (O reveals that the 
critical exponent s can be directly measured in an ex- 
periment by observing the dynamics of expansion of the 
atomic cloud at large times. The scaling theory of lo- 
calization predicts s = v for 3D disorder flil ]. A mea- 
surement of the critical exponents in an experiment with 
a Bosc-Einstein condensate could provide a spectacular 
test of the (one-parameter) scaling theory. Note that 
the decay of Sn with time predicted by Eq. (O is slower 
than what would have been obtained if localization were 
neglected (Sri ~ 1/t 3 ^ 2 ). Hence, localization effects mod- 
ify the dynamics of the condensate expansion even when 
e c < fJ- 

An important time scale can be obtained from the com- 
parison of Eqs. ([3]) and ([5]). Indeed, the convergence of 
the time-dependent profile fi(r, t) to fi(r, oo) given by Eq. 
© should become apparent when Eqs. ^ and ([5]) be- 
come of the same order. This happens at the "localiza- 
tion time" 



h fr\ 2 + s / v 



(0) 



If we accept that v = s in 3D [lfj, we obtain t\ oc oc r 3 
independent of v = s. Note that the localization time 
© exceeds the arrival time (@|, if r is larger than the 
healing length l/fc M of the initial condensate. 

We now consider the anomalously diffusing atoms with 
energies very close to the mobility edge. Indeed, in a nar- 
row energy strip | e — e c | < ecd^l/ec) 1 / 38 the critical diffu- 
sion, that led to Eq. ([5]), is taken over by the anomalous 



171 ]. Analysis shows that 



diffusion: D e (Q) oc (-iQ) 1 / 3 
for t > fioc the anomalously diffusing atoms dominate 
the dynamic part of the condensate density: 

-2/3-1/3S 



Sn(r,t) oc f( Cc /n)— (jt 



(7) 
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FIG. 1: Profiles of the average atomic density, associated with 
the localized part of the atomic cloud, at long times. The solid 
lines are obtained from the self-consistent theory of localiza- 
tion for three different ratios of mobility edge e c and chemical 
potential fi of the initial condensate. For each e c //x, the frac- 
tion of localized atoms /i oc is given. The dashed line shows 
the asymptote ft oc 1/r 4 followed by all curves for large r. 



However, at t ~ t\ oc the density of anomalously diffusing 
atoms (O is already of the same order as the density 
of localized atoms ([3]). Therefore, the contribution of 
anomalously diffusing atoms may be difficult to observe 
in an experiment. 

Up to here we have presented a general analysis based 
on scaling ideas only. The self-consistent theory of lo- 
calization [l|| allows us to calculate n(r, t) from Eq. 
(|2"|) without any additional assumptions. In this theory 
D € (Q) obeys 
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where g max = tt/3£. The self-consistent theory predicts 
the critical exponents v = s = 1. We plot n(r, oo) and 
Sn(r,t) obtained from Eqs. © and © in Figs. Q] and [2] 
As follows from Fig. [TJ for sufficiently large r the density 
profile n(r, oo) decays as a power-law 1/r 4 for all e c //i, as 
predicted by Eq. ([3]) with v = 1. The dynamic part of the 
atomic density Sfi(r, t) is shown in Fig. [2] and reaches a 
maximum at times of order of ^arrival given by Eq. ((4]) , at 
least when e c < fi. The maximum shifts to shorter times 
when e c //x increases. As predicted by Eq. ((5|), the plots 
of 5fi(r,t)r 3 /N for different r fall on a universal curve, 
when shown as functions of t/t aIT i va \ at fixed e c //i. All 
curves in Fig. [2] follow the 1/t asymptote for long times. 
This is in agreement with Eq. (O for t < t\ oc and Eq. ([7]) 
for t > t\ oc , with 8 = 1. 

The number of atoms that stay localized after the time 
t\ oc can be found by integrating Eq. ^ over r and by re- 
stricting the integration over energies to e < e c . We find 
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Einstein condensate expanding in a 3D random potential 
is governed by a single parameter equal to the ratio of 
the mobility edge e c to the chemical potential fi of the 
condensate. The critical exponents of the Anderson lo- 
calization transition v and s determine the evolution of 
ri(r, t) in time and space. Our results open a new way to 
measure the critical exponents in an experiment. 
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FIG. 2: Profiles of the dynamic part of the average atomic 
density obtained from the self-consistent theory of localiza- 
tion (solid lines). The time is given in units of tarrivai de- 
fined by Eq. (jlj). The dashed lines show 1/t asymptotes. At 
a given ratio of mobility edge e c and chemical potential fj,, 
5n(r, t)r 3 /N for different r (different colors) fall on a univer- 
sal curve. For clarity, the curves corresponding to e c /fi = 1 
and 10 are shifted downwards by 1 and 2 units, respectively. 
The inset shows the complete atomic densities n(r, t) obtained 
by adding Sn shown in the main plot and n(r, oo) of Fig. [T] 



for the fraction of localized atoms /i oc ~ (5/2)(e c //i) 3//2 
for e c <C [i- and f\ oc ~ 0.45 for e c \x. Remarkably, 
this result is independent of the model used for the lo- 
calization length £(e) and, in particular, independent of 
the critical exponent v. The conclusion that even for 
strong disorder (e c ^> fx) only a fraction of the conden- 
sate stays localized is a surprising outcome. It is due 
to the large-energy tail of the spectral function derived 
above for the uncorrelated disorder: A(k, e) cx 1/e 3 / 2 . 
Atoms with high energies exhibit diffuse behavior, even 
though their initial kinetic energies were small (ek <C /Lt). 
This explains the small value of /i oc . It is straightfor- 
ward to show that for a potential with correlation radius 
ro the large-energy tail of A(k, e) is suppressed for ener- 
gies larger than eo ~ e c (£/r ) 2 . Thus, correlations seem 
to help in localizing more atoms. 

The available ID experiments Q and theoretical devel- 
opments [TH allow us to estimate the minimal realistic 
values of the mean free path I in 3D optical speckle po- 
tentials to be in the range of 1 to 10 fim. From this 
values we obtain e c /fi ~ 10~ 4 -10~ 2 for 87 Rb atoms and 
fj,/h ~ 10 kHz @. The corresponding arrival (localiza- 
tion) time is in the range of 1 to 10 ms (100 ms to 100 s) 
for r = 5£. With the lifetimes of Bose-Einstein conden- 
sates attaining 20 s 



2l|, these estimations demonstrate 



that our theoretical predictions are within the reach of 
current experiments. 

In conclusion, we have shown that the long-time large- 
distance behavior of the average density n(r, t) of a Bose- 
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